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We investigate single-inclusive high-pr jet production in longitudinally polarized pp collisions at 
RHIC, with particular focus on the algorithm adopted to define the jets. Following and extending 
earlier work in the literature, we treat the jets in the approximation that they are rather narrow, in 
which case analytical results for the corresponding next-to-leading order partonic cross sections can 
be obtained. This approximation is demonstrated to be very accurate for practically all relevant 
situations, even at Tevatron and LHC energies. We confront results for cross sections and spin- 
asymmetries based on using cone- and fct-type jet algorithms. We find that jet cross sections at 
RHIC can differ significantly depending on the algorithm chosen, but that the spin asymmetries are 
rather robust. Our results are also useful for matching threshold-resummed calculations of jet cross 
sections to fixed-order ones. 
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J2' 



> 

in 
00 

OV 

o: 

(N 



X 



I. INTRODUCTION 

Jets are copiously produced at high-energy hadron col- 
liders. Among other things, they play important roles as 
precision probes of QCD and nucleon structure. At the 
Relativistic Heavy Ion Collider (RHIC), jets are by now a 
well-proven tool for investigating the spin structure of the 
nucleon, through double-helicity asymmetries measured 
in the reaction pp — > jet X. The corresponding measure- 
ments [l[ have in particular provided exciting information 
on the proton's polarized gluon distribution, Ag. 

There is no unique way of defining a jet. As a re- 
sult, different jet algorithms exist and are being used in 
experiment. The jet definitions and algorithms can be 
broadly divided into two classes [3] : (i) successive recom- 
bination algorithms |h3, and (ii) cone algorithms 0. 
For the former, one first defines a distance between a 
pair of produced objects (initially, two particles), as well 
as a "beam distance" of each object to the collider beam 
axis. For each object, the smallest of these distances is 
determined. If it is a beam distance, the object is called 
a jet and removed from the list of objects in the event; 
otherwise the two objects are combined into a single ob- 
ject. This procedure is repeated until no further recom- 
binations take place. Prominent examples of successive 
recombination algorithms are the so-called k t [1,|I[, Cam- 
bridge/Aachen 0], and anti-fct @ algorithms, which dif- 
fer in how the distances are defined. We will collectively 
refer to such algorithms as "fet-type" algorithms. 

Cone algorithms also come in different variants 
They have in common that the jet is defined by the par- 
ticles found inside a circle in the plane formed by ra- 
pidity and azimuthal angle, such that the sum of the 
four-momenta of these particles points in the direction of 
its center. While widely used in experiment, the tradi- 
tional cone algorithms (notably the ones known as "mid- 
point" algorithm mostly used at RHIC [l|, Q and the "it- 
erative cone" algorithm) were found to be not infrared- 



safe [lfl, [ll|. This evidently sets a serious limitation 
to the use of such algorithms in the theoretical calcula- 
tion and to comparisons of data and theory. For single- 
inclusive jet cross sections, the lack of infrared-safety 
becomes an issue first at next-to-next-to-leading order 
(NNLO) in perturbation theory, so that next-to-leading 
order (NLO) calculations remain meaningful in the sense 
that they produce finite and well-defined answers. In 
case of the midpoint cone algorithm, a solution to the 
problem of infrared-unsafety was found in terms of the 
"Seedless Infrared Safe Cone" (SISCone) algorithm [111 ]. 
It was also shown that the anti-fc t algorithm mentioned 
above can effectively cure the lack of infrared-safety of 
the iterative cone algorithm [fj. As a result, the SIS- 
Cone and all fc t -type algorithms are nowadays known to 
be infrared-safe and are preferred for use in experiments. 

In earlier work [13, [H| , the spin-dependent (and spin- 
averaged) cross sections for pp — > jet X were derived at 
NLO. Ref. [13] was based on a Monte-Carlo integration 
approach, while 13] used a largely analytic technique for 
deriving the relevant partonic cross sections for cone al- 
gorithms, which becomes possible if one assumes the jet 
to be a rather narrow object [r3T4l7l ]. This assumption 
is equivalent to the approximation that the cone open- 
ing R of the jet is not too large, and hence was termed 
"Small Cone Approximation" (SCA) in [13]. In the SCA, 
one systematically expands the partonic cross sections 
around R = 0. The dependence on R is of the form 
A\ogR + B + 0(R 2 ). The coefficients A and B are re- 
tained and calculated analytically, whereas the remain- 
ing terms 0(R 2 ) and beyond are neglected. The advan- 
tage of the analytical method is that it leads to much 
faster and more efficient computer codes and is hence 
readily suited for inclusion of jet spin asymmetry data 
from RHIC in a NLO global analysis of polarized parton 
distributions. Indeed, the results of [l3[ have been used 
in the global analysis fisl ] . where the experimental data 
from STAR [l| were used to constrain Ag. It was shown 
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in |l3[ that the SCA is in fact an excellent approxima- 
tion to the full Monte-Carlo calculation (which is valid 
for arbitrary cone openings), for all values of R and kine- 
matics relevant at RHIC. (We shall revisit this finding in 
our phenomenological section Hill below) . 

Commensurate with the procedure chosen by STAR, 
the calculation [l3] was performed for the midpoint cone 
algorithm. In the present paper, we will extend the work 
in [l3[ to the case of fct-type algorithms. We will again 
use the approximation of a rather narrow jet. As we do 
not really have a jet "cone" for the fct-type algorithms, we 
shall from now on refer to this approximation as "Narrow 
Jet Approximation" (NJA) . This term will be collectively 
applied to both the cone (where it used to be the SCA) 
and the fct-type algorithms. The meaning of the "NJA" 
will always be that the jet parameter R used to define 
the jet (cone opening for the cone algorithm and distance 
between two objects for the fct-type algorithms) is not too 
large, as we shall discuss in more detail below. 

One motivation for our new study is that fct-type al- 
gorithms are also being considered in STAR now (l9j . 
so that it is timely to prepare the corresponding theo- 
retical NLO calculations for the spin asymmetries. The 
differences between the jet cross sections for the cone and 
fct-type algorithms in the NJA are also interesting from a 
theoretical point of view. We will find that they amount 
to finite contributions with leading order (LO) kinemat- 
ics. As such they play a role as matching coefficients in 
threshold resummation studies of jet production, as was 
discussed in (2^] . They also appear in a related context in 
studies of jet shapes in the framework of "Soft Collinear 
Effective Theories" (SCET) jp. 

The remainder of this paper is organized as follows: 
In Sec. HH we present the technical details and analyti- 
cal results of our calculation of single-inclusive jet cross 
sections in the NJA, focusing on the fct-type algorithms. 
Section lTOl contains phenomenological results relevant for 
RHIC. We summarize our work in Sec. IIVI 



For the fct-type algorithms one defines for each 

pair of objects (initially, particles) j, fc the quantity 



Ijk 



r>2 



(2) 



where p is a parameter that specifies the algorithm, kx 
denotes the transverse momentum of particle j with re- 
spect to the beam direction, and 



(Vj ~ Vk) 2 + {fa ~ 0fc) 2 



(3) 



The parameter R is called the jet radius, djk may be 
viewed as a distance between two objects j and k. One 
also defines for each object a distance to the initial beams: 



IjB 



k 2p 



(4) 



The algorithm identifies the smallest of the djk and djB- 
If it is a beam distance, the object is defined as a jet and 
removed from the list of objects. If the smallest distance 
is a djk, the two objects j, k are merged into a single 
one. The procedure is repeated until no objects are left 
in the event. As mentioned above, the jet algorithm is 
fully specified by the parameter p. We have p = 1 for 
the fct-algorithm [3,0], P = for the Cambridge/ Aachen 
algorithm 5], and p = — 1 for the anti-fct algorithm |6[. 

Note that on top of the choice of jet algorithm one also 
has to define how objects are to be merged if the need for 
that arises. Throughout this paper our choice will be (for 
both algorithms) to define the four-momentum of a new 
object as the sum of four-momenta of the partons that 
form the new object. This so-called "E recombination 
scheme" 0] is the most popular choice nowadays. 



B. Calculation of single-inclusive jet production 
cross sections at NLO 



II. TECHNICAL DETAILS 



The spin-averaged cross section for the process 
p(P a )p(Pb) — > jet(pj) X can be written as [13| 



A. Cone and fct-type jet definitions 

We consider single-inclusive jet production in hadronic 
collisions, pp — > jetX, where the jet has a transverse mo- 
mentum ptj, pseudorapidity rjj, and azimuthal angle <pj. 
The cross section is infinite unless a finite jet size is im- 
posed as a parameter. The different jet algorithms vary 
in the way this size is defined. In the cone algorithm [7[ , 
one defines the jet by all particles j that satisfy 



(vj-v,) 2 + (<l>j-<Pi) 2 <R 2 - 



(1) 



Here rjj and cj)j denote the pseudo-rapidities and az- 
imuthal angles of the particle, and R is the jet cone 
aperture. The jet four-momentum sets the center of the 
cone; it is nowadays usually defined as the sum of the 
four-momenta of the particles j forming the jet. 



E 



dv 



dpTjdrjj S J vw v(l - v) 



dw 



< I fa(x ai fX F )f b (x b ,^ F ) 

'VW/v w 



da 



(o) 

ab— ^jetX 



(s,v) 



dv 



5(1 -w) + 



dvdw 



,(5) 



where the dimensionless variables V and W are defined 
in terms of pxj and rjj as 



PTj 



and W- 



(6) 
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with y/S = \/ (P a + Pij) 2 the hadronic c.m.s. energy, v 
and w are the corresponding parton-level variables; they 
are given in terms of the partonic Mandelstam variables 

S=(p a +Pb) 2 , t= (Pa-Pj) 2 , U=(p b -pj) 2 , (7) 



as 



i + i, 

s 



s + t 



(8) 



The sum in ([5]) runs over all partonic channels a 



do) 



and da 



the LO and NLO 



jet + X, with da ab ^ ctx ^ ^ ab ^ ]c} x 
terms in the corresponding partonic cross sections, re- 
spectively. f a (x ai Hf) and fb(xb,p.F) denote the par- 
ton distribution functions at factorization scale p,F whose 
partonic momentum fractions are determined by V, W, v 
and w: 



VW 



1 - V 
1-v' 



(9) 



Finally, [Ar in (JSJ is the renormalization scale for the 
strong coupling constant. 

We note that expression (J5J can be straightforwardly 
extended to the case of collisions of longitudinally po- 
larized protons. Here one defines a spin-dependent cross 
section as 



d 2 Aa 1 
dpTjdrjj ~ 2 



d 2 a++ d 2 a^ 



dpTjd-qj dpTjdrjj 



where the superscripts indicate the helicities of the col- 
liding protons. The structure of Eq. ([5]) also applies 
to d 2 Aa/dpT,dr]j, except that the partonic cross sec- 
tions and parton distribution functions are to be replaced 
by their spin-dependent counterparts dAa a b^,j C tx and 
Af a ,b(xa,bi Hf), respectively. The former are defined in 
analogy with (TXOl) . and for the latter we have 



(11) 



Af a (x,fj, F ) = f+(x,/j, F ) - f a (x,ix F ), 



where /+ (/~) denotes the distribution for partons of 
type a with the same (opposite) helicity as that of the 
parent proton. All our expressions below will be formu- 
lated for the spin-averaged case; however, they equally 
apply to the polarized case with the modifications just 
discussed. 

A possible way of organizing the NLO calculation of 
the single-inclusive jet cross section was developed and 
employed in Refs. [H, [TM3, El- It starts from the 
NLO single-parton inclusive cross sections da a b^ c x , rel- 
evant for the single-inclusive hadron production process 
pp — > hX and analytically known from previous calcu- 
lations [H, [23[ . These cross sections cannot directly be 
used to describe jet production; they can, however, be 
converted to the desired single-inclusive jet cross sections. 
To this end, one first imagines a "jet cone" around the ob- 
served parton c and notices that a NLO single-parton in- 
clusive cross section contains configurations where there 
is an additional parton d inside the cone (note that we 



use the term "cone" here just for simplicity- the con- 
siderations apply to any jet definition). For a jet cross 
section, the observed final state should not just be given 
by parton c, but by partons c and d jointly. One therefore 
subtracts these contributions and replaces them by terms 
for which partons c and d are both inside the cone and 
form the observed jet together. To be more precise, for 
a given partonic process ab — > cde we have, after proper 
bookkeeping of all partonic configurations that are pos- 
sible in the cone: 



d<Jab^jctX = [d<T c — d& c ( d ) — rf(7 c ( e ) 

+ [da d - rf<7 d ( c ) - da d ( e } 

+ [da e 

+ d& cd 



do e {c) - da e (d)} 
da ce + dade- 



(12) 



Here d&j is the single-parton inclusive cross section where 
parton j is observed (which also includes the virtual cor- 
rections), d&j^ is the cross section where parton j is 
observed but parton k is also in the cone, and dajk is the 
cross section when both partons j and k are inside the 
cone and jointly form the jet. 

The single-parton inclusive cross sections d&j of [l6l . l23l ] 
were obtained after a subtraction of final-state collinear 
singularities in the modified minimal subtraction (MS) 
scheme. Upon calculation of the combinations daj^ + 
(10) dak(j) ~ ddjk in (jf 2[) one also finds collinear singularities, 
which must match those initially present in (d&j +dak)/2- 
On the other hand, the full expression in Eq. (fT2")l . be- 
ing an inclusive-jet cross section, must be collinear-finitc. 
Therefore, in order to obtain the combination in (fT2"]l cor- 
rectly, one just needs to perform an MS subtraction also 



da 



Ho) 



da 



For 



the situation is different since here bo th p articles j and 
k jointly form the jet. In Refs. [TJ, UH, Il7l | the d&jk were 



of the singularities in the d&j^ 
further discussion, see [l3j . 

In practice, it is convenient to consider the d&j^) and 
d&jk separately. In the NJA, to which we will turn in the 
next subsection, they may in fact be computed analyti- 
cally. At NLO, they both receive contributions only from 
real-emission 2 — > 3 diagrams. Since for d&j^) the jet is 
obtained from the single parton j , it is independent of the 
jet definition. The relevant results for this piece in the 
NJA may be found in [f3j]. For d&jk, on the other hand, 

here bo th p a 
[II El El t 

obtained for the case of cone algorithms. In the following 
we will derive them also for the fc t -type algorithms, using 
again the NJA. 



C. Calculation of dajk in the "Narrow Jet 
Approximation" 

It is instructive to discuss the cone and fct-type al- 
gorithms in parallel, in order to make contact with the 
derivations made in [f3[ and to make our paper self- 
contained. From Sec. Ill Al we see that both the cone and 
the fct-type algorithms contain a jet parameter R. In the 
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NJA one assumes R to be relatively small. As discussed 
at the end of Sec. Ill Al our choice is to merge objects 
by adding their four-momenta. For dajk this means that 
the four-momentum pj of the jet is the sum of the parton 
four-momenta pj and pk ■ 

We first observe that for the fej-type algorithms the 
two partons j, k are merged into one jet if their distance 
defined in @ is smaller than their respective beam dis- 
tances diB and djB defined in (j4|). For dajk this has to 
be the case by definition, and we therefore arrive at the 
condition 



for fet— type algorithms, 



(13) 



with Rij defined in Eq. We stress that this condition 
holds for all fc 4 -type algorithms, regardless of the choice 
of the parameter p. This implies that at NLO all fcj-type 
algorithms lead to the same jet cross section, a result 
that does not rely on the NJA. Equation (TT3"f is to be 
contrasted with Eq. (TTJ) valid for the cone algorithm: 



R]j < R 2 A Rlj < R 2 



for cone algorithm. (14) 



The difference between Eqs. (IT31 and (TTJ} is solely re- 
sponsible for any differences between the NLO results 
for the two types of jet algorithms! For the fet- type algo- 
rithms it is the "distance" between the two partons that 
is constrained by the jet algorithm, whereas for the cone 
algorithm it is the distance of each parton to the jet itself. 
Note that this observation was already made in Ref. |3j, 
and in 21] in a slightlydifferent context. 

As was shown in [ljj], in the NJA dajk is, up to 
trivial factors, given by the following expression (see 
Eqs. (19), (20), (27) of that paper): 



da 



dvdw 



dPS 3 Pf K (z) 
dvdw 2pj ■ pk ' 



(15) 



where the integration dPS 3 is over the phase space of 
the three-body final state of the overall partonic pro- 
cess, which is carried out in d = 4 — 2e dimensions. Ex- 
pression (|15p arises from the fact that dajk is strongly 
dominated by contributions for which particles j and k 
result from collinear splitting of an intermediate parti- 
cle K. The reason is that in this case the propagator of 
the intermediate particle, represented by the denomina- 
tor l/(2pj ■ Pk), goes on-shell. For instance, if the jet is 
formed by a quark and a gluon, the pair will predomi- 
nantly originate from a quark splitting into a quark plus 
a gluon, described by the splitting functions P qq and P gq . 
The argument z of the splitting function is the fraction 
of the intermediate particle's momentum transferred in 
the splitting. The superscript "<" on the splitting func- 
tion indicates that the d-dimensional splitting function 
Pjk{z) is strictly at z < 1, that is, without its 6(1 — z) 
contribution that is present when j = K . This is a nec- 
essary condition in order to have two partons producing 
the jet. For additional details, we refer the reader to 
Ref. [H. 



Making use of the fact that pj = pj + Pk, the term on 
the right-hand-side of dT5|) may be written as fl3T | 



dPS 3 Ff K {z) 




d0^ 



l-2e , 



COS 9jk 



(16) 



where Ej = Ej + Ek is the jet energy (with Ej k the 
energies of partons j, k), z — Ej/Ej, and 6jk the angle 
between the three- momenta of partons j and k. Qj is the 
polar angle of parton j, measured with respect to the jet 
direction. # max is an upper limit on the -integration 
that needs to be derived according to the jet algorithm. 
It is of the order of the jet parameter R, and hence treated 
as small in the NJA. 

It is useful to write the ^-integral as an integral over 



the (squared) invariant mass p 2 j 
produced jet. One finds 



2pj ■ pk of the 



s(0 jk) 



cos 



/) = 



2EjEk 
2EjEj - m 2 
2E^/Ej^?, 



(17) 



With this one obtains after some algebra: 



dPS 3 Pf K (z) 
dvdw [2pj ■ pk) 



47ry e (v(l-v))- 
64tt 3 \y/s) r 2 (l-e) 



xS(l-w) / dzz- £ (l - z)- e P< K {z) 



dm 2 



-2e 



(18) 

where we have also expressed the integral over the energy 
Ej as an integral over z. In the NJA, m 2 is a small quan- 
tity, and we have hence neglected powers of m 2 wherever 
possible. Note however that the m 2 -integral produces a 
1/e-singularity at the lower end. 

All that is left to be done now is to determine the 
upper limit of the integral over m 2 , which depends on the 
jet algorithm chosen. In order to make contact with the 
results in Ref. [IH , we do this first for the cone algorithms 
and afterwards for the fc t -type algorithms we are mainly 
interested in here. In both cases we first write the jet 
four-momentum as 



Pj 



where 



Ej, \pj 



r) 



\P 



cosh(?7j) cosh(?7j) 



, |p i7 |tanh(77 J 



\pJ 



E 2 



(19) 



(20) 
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For the parton momenta pj and pk, which are light-like, 
we write accordingly 

_ / cos(6j) sin((/),,) , , A 

« = £ H 1 '^'^' tanh H' 



'101 "vv-VHs 

cos(<fa) sin(0 fc ) 
cosh(ry fc ) ' cosh(?7 fc ) 



,tanh(77 fc ) , (21) 



with the azimuthal angles (f>j^ and pseudorapidities T)j t k 
of the partons. 

(i) Gone algorithms: We write 

to 2 = 2pj ■ p k = 2pj ■ pj 

cosh(?7 i - rjj) - cos(0j - (j>,j) 



2E,E 



cosh(7jj) cosh(?7j) 



Ejm 2 cos(0j — </>j) + sinh(ryj) sinh(?7j) 



cosh(77j) cosh(?7j) 



(22) 



where we have expanded \pj\ to first order in m 2 . The 
combination (cosh(v7j — 77,7) — cos(0j — 0,/)) in (j22j) is 
small, while in the other term m 2 is small. Therefore, 
to the order we consider, we can set rjj — nj, <pj = tpj in 
all other places. This gives 



E,,E. 



cosh (77 j ) 



Ej 



(23) 



where Rjj is as defined in Eq. ((T|). We solve for m 2 and 
get 



cosh (?7j) 1-2 



Likewise, we find 



1 - z 



cosh (77,7) ^ 



1? 



fej' 



(24) 



(25) 



The jet criterion (1141) in the cone algorithm then imme- 
diately translates into 



E 2 R 2 



max, cone 



l-z 



(26) 



cosh z (?7j) V 1 ^ 2 2 
The last two integrals in (|18[) are now readily performed: 

dm 2 



dzz- £ {l- z)- £ P<{z) 



-2e 



1 / ff 2 i? 2 



£ \cosh (77 j) 



(27) 



where 



l jK 



1/2 



dzz 



-2e 



+ / dz(l - z)~ 
1/2 



(28) 



The explicit results for these i nteg rals for the various 
splitting functions were given in [131 ]: 



rconc 
99 



rconc 
99 



rconc 
ff9 ' 



1 



7 7T 2 

2 + T 



3 log 2 



1 

2 

2C. 



2 /23 
3 + Hl8 



log 2 



1 11 



137 7T 2 11, „ 

^6- + y-y log2 



(29) 

with Ca — 3 and Cf = 4/3 the usual SU(3) Casimir 

operators. Note that the ratio R/ cosh(?7j) corresponds 

to the half-opening S of the jet cone considered in [l3j]. 

We also note the logarithmic dependence of d&jk on R 

arising from the factor R~ 2e in Eq. ([27]). 

(zi/ 1 kt-type algorithms: Here we use Eq. (|2ip to directly 

compute 



2EjEk 



(cosh(77 3 - - %) - cos(<£,- - ^> fe )) 



cosh(7yj ) cosh(r/fc) 
We can approximate this by 

((% - %) 2 + " ^) 2 ) 



(30) 



EjEk 



cosh (77,7) 
EjEk 



cosh 2 (77 j 



■I? 



(31) 



with i? jfc defined in ([15]). The condition i? 2 fc < R 2 then 
immediately gives 



777 



cosh (77 j ) 
and instead of (|2"T)) we have 

el 

d««- e (l-«)- e P^(z) 



(32) 



• fc t dm 2 



,-2e 



1 / E 2 R 2 
e V cosh 2 (77 j ) 



jk t 



(33) 



where 



^9 = Cf 



£ - 



13 2tt 2 

T + — 



rfet 
99' 



jfc* 
?3 



T 

99 



2C A 



23 



1 11 

"e ~ ~6~ 



67 27T 2 ; 
9 + 3 



(34) 



6 



Comparison with (|29l) shows that the pole terms in 
Eqs. ([2"?} . (j3"3")l are the same, as they have to be. The 
finite remainders, however, differ and will lead to finite 
and i?-independent differences in the NLO cross sections 
for the two types of algorithms. As seen from Eq. (|T6"|) . 
the cross sections dajk are proportional to 6(1 — w) and 
hence have LO kinematics. This will, therefore, also be 
true for the finite differences just mentioned. We note 
that expressions similar to those in (j3~4")l were also ob- 
tained in [2l| in the context of jet studies in SCET. 

It is now straightforward to use the integrals I^k given 
above to derive the NLO jet cross section for the fct-type 
algorithms in the NJA- the calculation proceeds exactly 
as in |13lLln fact, is is very easy to change the numerical 
code of [13| to the case of the fcj-type algorithms: The 
integrals Ij]^ e are the only sources of terms oc log 2 in 
the NLO calculation for the cone algorithms. Thus, by 
replacing these terms appropriately in each subprocess 
one can translate the result from the cone algorithm to 
the case of fct-type algorithms without having to do the 
full calculation in (H"2"1) . 

While we have derived all results above for the spin- 
averaged cross section, it is straightforward to extend 
them to the case of jet production in polarized collisions. 
In the NJA, the contributions by particles j and k form- 
ing the jet entirely arise from final-state emission, which 
is independent of the polarization of the initial partons. 
Therefore the same integrals ij™ c or l-^ apply to the 
polarized case. 



III. PHENOMENOLOGICAL RESULTS 

In this section, we present some phenomenological re- 
sults for single-inclusive jet production cross sections and 
spin asymmetries in pp collisions at RHIC. Our main fo- 
cus is of course on the sensitivity of these quantities to 
the jet algorithm used. 



A. Unpolarized collisions 

We begin by ascertaining the accuracy of the NJA. It 
was shown already in [l| that for the cone algorithm the 
approximation is very accurate for the values of R and 
kinematics relevant at RHIC. To confirm this finding, we 
make use of the recently developed fastNLO package [24j 
which is based on the NLO code NLOJet++ of [25[ and 
allows to compute NLO jet cross sections with Monte- 
Carlo integration methods. In fact, the authors of the 
code offer an online tool that provides numbers for the jet 
cross section at RHIC for the mid-point cone algorithm, 
as used by STAR In the upper left part of Fig. [T] we 
compare these results to the ones we find with our code 
based on the NJA. We plot the ratio of the two theoret- 
ical results. We have used here the CTEQ6.6M parton 
distributions [26| . which will be our choice for the spin- 
averaged parton distribution functions throughout this 



paper. The comparison is for y/S = 200 GeV, R = 0.4, 
and a range of pseudorapidity 0.2 < |?7j| < 0.8. We 
have chosen the factorization and renormalization scales 
as [if = Mi? = Pt,j- As one can see, there is excellent 
agreement between the full NLO Monte-Carlo calculation 
and our approximated result. The largest deviations oc- 
cur at the lowest ptj ; even here they amount to at most 
5%. 

The lower left part of Fig. [T] shows a similar compari- 
son for the fct-type algorithms (we remind the reader that 
at NLO the jet cross section is the same for all fcj-type 
algorithms). The exact NLO calculation was performed 
here with the Fast Jet code [13, [H[ , which is also based 
on (25[. Kinematics are similar as before, except that we 
have used here the rapidity range |t7j| < 0.6. We have 
again used R = 0.4 for the jet parameter in the fct-type 
algorithms. Again, the NJA reproduces the full NLO 
calculation very accurately. Interestingly, comparing the 
upper and lower left parts of Fig. [TJ we observe that 
the NJA very slightly overpredicts the NLO cross sec- 
tion for the case of cone algorithms, but underpredicts it 
in the fc t -type case. We note that the excellent overall 
agreement between the NJA and the exact NLO calcula- 
tion occurs also for other choices of the factorization and 
renormalization scales. 

As an aside, we also show in the right part of Fig.[T]re- 
sults for similar comparisons for pp collisions at the Teva- 
tron (VS = 1960 GeV) and for the LHC (y/S = 7 TeV), 
for various jet algorithms. The kinematics correspond 
to those used in experiments [29|-[32j]. The exact NLO 
results wer e ag ain obtained using the fastNLO [24j and 
FastJet [27, 28] packages. One can see that the NJA also 
works very well in these cases. 

Another way of gauging the accuracy of the NJA is 
to consider the ratio of jet cross sections for different jet 
parameters R: 



n{R 1 ,R 2 ) 



d 2 a/dpTjdrij(R = Ri) 
d 2 a/dpTjdi]j(R = R 2 ) 



(35) 



As shown in [33j |. TZ can be expanded perturbatively in 
orders of a s . To the lowest non-trivial order one has 

^^NLO/n \ j2 NLO ( n \ 

■uiu u\ 1 i a (Ri)-aa (R 2 ) , 
n(R 1 ,R 2 ) = 1 + „ att q I , (36) 



rf2(T NLO 



where d 2 cr NLO (i?) denotes the NLO cross section for a 
given R and d 2 a Nh ° \o(a 2 ) its truncation to the lowest 
order, keeping however the strong coupling constant a s 
and the parton distributions at NLO. d 2 a Nh \o(a 2 ) does 
not depend on R. The difference of cross sections in the 
numerator of (|31)|) is of order a^, so that 1Z(R\,R 2 ) is 
of the form 1 + 0(a s ). Figure [2] shows the result for 
ft(0.2, 0.4) at RHIC energy \/S = 200 GeV, as a func- 
tion of pt_, ■ The cross sections have been integrated over 
\vj\ < 0.6, and we have used here scales [ip — — Pt.j- 
Our result may be directly compared to the correspond- 
ing one given in (3^1 for the same set of parameters, also 
shown in the figure, where the full FastJet code was used. 
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One can see that the agreement is excellent, impressively 
demonstrating the accuracy of the NJA. We note, how- 
ever, that in the NJA the ratio 7Z(Ri,R2) is independent 
of the jet algorithm chosen. The agreement seen in Fig. [2] 
hence is a test of the NJA as such, but not of the imple- 
mentation of a specific jet algorithm. 

Having established the validity of the NJA, we now 
provide results for jet cross sections at RHIC. Figure [3] 
shows the spin-averaged cross sections for |?7j| < 1 at 
VS = 200 GeV (left) and y/S = 500 GeV (right). Re- 
sults are presented for both the cone and the fc t -type al- 
gorithms, using two values for the jet parameter, R = 0.4 
and R = 0.7. The renormalization and factorization 



scales have again been set to pxj ■ Figure 0] examines 
how the cross sections vary with the choice for the scale 
[i = iip = fix, for the case R — 0.4. We vary the scales 
in the region pt. 7 /2 < < 2pTj and show the relative 
deviation from the result for the case /i = pTj • Interest- 
ingly, for this value of R, the scale dependence is not too 
similar for the cone and the fct-type algorithms. For the 
former, the cross section increases monotonically when 
going from scale [i = 2pxj to \i = Pt,/2, while for the 
latter the cross section for /x = px a is for most px, larger 
than those for both other scales. Also, the scale depen- 
dence is overall somewhat smaller for the fct-type algo- 
rithms. We have verified that these patterns are present 
in the exact NLO calculation with FastJet [2o| . 

We now turn to a more detailed comparison of the 
jet cross sections at RHIC for the two different jet algo- 
rithms. We define the ratio 



[d 2 (A)a/d PTj drj 



J i fc t -type 



[d 2 (A)a/dpTjdr]j] 



(37) 



choosing the same jet parameter R for both cross sec- 
tions. Figure [S] shows our results for TZ a \ so as a function 
of pt, (we have chosen pr, bins of 5 GeV width), for 
R = 0.4 and 0.7. We present results for both energies rel- 
evant at RHIC, VS = 200 GeV (left) and VS = 500 GeV 
(right). We have in both cases integrated the cross sec- 
tions the over the pseudorapidity range 1 77^7 1 < 1. We 
have computed the ratio for the scale \x — PT. n keeping 
in mind, however, that according to Fig. 2] its scale de- 
pendence is quite large. As one can see, the cross section 
for the fct-type algorithms is about 10% smaller than that 
for the cone algorithm, except at pxj ^ 10 GeV where 
the ratio 7?. a i go drops more strongly. Our results are con- 
sistent with the trend seen in jet algorithm studies by 
STAR 9]. The ratio 7\L a ig also shows relatively little 
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dependence on the jet parameter R. 

The fact that the jet cross section for the fct-type al- 
gorithms is found to be smaller than that for the cone 
algorithm for the same value of R implies that a choice 
of a larger R for the fc(-type algorithms should bring the 
two cross sections much closer together. Indeed, it was 
found in Ref. Q, that a choice Rk t ~ 1.35i? con c makes 
the cross sections for the two algorithms quite similar, 
also for other choices of the scale /i. We confirm this 
finding. 

B. Longitudinally polarized collisions 

For the polarized case we use the "DSSV" helicity par- 
ton distributions of Ref. [18]. Our first finding is that for 
the polarized case the effects of changing the jet algo- 
rithm are somewhat more pronounced than in the unpo- 
larized one. Figure [S] shows the ratio 7?. a igo for polarized 
collisions at RHIC, again computed for the scale fi — pxj ■ 
72-algo is again around 90% at high pT r , but shows large 
deviations from unity in the bin around pt, = 12.5 GeV. 



The reason for this is that for the DSSV set of parton 
distributions the polarized jet cross section crosses zero 
around pt, = 10 GeV. Depending on the jet algorithm, 
this zero will be at slightly different locations, making 
the denominator and numerator of 7?. a igo vastly different 
there. This is, of course, for the most part an artifact 
of the way we are performing the comparison of the jet 
cross sections, taking ratios of small numbers at some 
Pt, ~ 10 GeV. On the other hand, it does demonstrate 
the issue that in regions where the polarized cross section 
is very small it may also be quite susceptible to the choice 
of jet algorithm and hence (at the non-perturbative level) 
to hadronization corrections. We note that the difference 
between the cross sections for the fct-type and cone al- 
gorithms may again be diminished by choosing a larger 
value of R for the former, R^ t ~ 1.35-R cone a s in the 
spin- averaged case discussed above. This also brings the 
two polarized cross sections somewhat closer together in 
the bins near their zero, even though marked differences 
remain here. 

Figure [7] shows the spin asymmetries All at RHIC, 
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which are defined by 



.4 



LL 



_ d 2 Aa/dpTjdrjj 
d 2 a/dp Tj dr]j 



(38) 



For the denominator we use the spin-averaged cross sec- 
tions shown in Fig. [3] The most important observation 
is that the asymmetries arc quite insensitive to the jet 
algorithm chosen, and also to the value of the jet param- 
eter R. The exceptions are of course regions where the 
polarized cross section (nearly) vanishes, as we saw in 
Fig. [5J In these regions, All is very small, and so these 
exceptions are not really noticable in Fig. [7J 



IV. CONCLUSIONS 

We have computed the NLO cross sections for single- 
inclusive high-pT jet production in spin-averaged and lon- 
gitudinally polarized pp collisions at RHIC, with special 



focus on the algorithm adopted to define the jets. Follow- 
ing Ref. [l3| , we have treated the jets in the approxima- 
tion that they are rather narrow ( "Narrow Jet Approxi- 
mation"). In this approximation one can derive analyti- 
cal results for the corresponding partonic cross sections, 
which are of the form A log R + B with R the jet pa- 
rameter. We have extended the results of Ref. [13[ to the 
case where an infrared-safe "A; t -type" algorithm (fc t , anti- 
kt, Cambridge/ Aachen algorithm) is used. By compari- 
son to available "exact" NLO jet codes for spin-averaged 
scattering [24], [27j , we have found that the Narrow Jet 
Approximation is very accurate at RHIC for practically 
all relevant situations. The same is true even at Tevatron 
and LHC energies. 

Our numerical results show that, for given R, jet cross 
sections at RHIC depend significantly on the algorithm 
chosen. Moreover, the scale dependence of the cross sec- 
tions can be quite different for cone- and fct-type algo- 
rithms. For polarized cross sections, the dependence on 
the jet algorithm can be very pronounced in the vicinity 
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of jet cross sections to fixed-order ones. For the case of 
cone algorithms, based on the results of [lj|, this was 
already exploited in [20j. Our present calculation allows 
to extend this procedure to the case of the nowadays 
more popular A: t -type algorithms. We note that the jets 
we consider here remain massive near partonic threshold 
(see the discussion in which affects the logarithmic 

structure of the partonic cross sections [34| and corre- 
sponds to the situation encountered in experiment and 
in the "exact" NLO codes such as FastNLO and Fast- 
Jet. It is known that "non-global" logarithms arise in 
this case 1351] . 



P T (GeV) 



FIG. 7: Double-longitudinal spin asymmetries All at RHIC, for 
y/S = 200 GeV and -\/S = 500 GeV and various jet definitions. 
We have averaged over \r]j\ < 1. The scales have been chosen 

as flp = fJ-R = PTj ■ 

of a zero of the cross section. On the other hand, spin 
asymmetries at RHIC overall turn out to be quite robust 
with respect to the jet algorithm adopted. 

We finally stress that our analytical results are also 
relevant for matching threshold-resummed calculations 
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